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ABSTRACT 

We use the semi-analytical approach to analyze gravitational lensing of remote 
quasars by foreground dark halos in various cold dark matter (CDM) cosmolo- 
■ gies, in order to determine the sensitivity of the predictions for probabilities of 

images separations to the input assumptions regarding the properties of halos 
. and cosmological models. The power spectrum of primordial fluctuations is nor- 
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malized by the cluster abundance constraints. The mass function of dark halos 



is assumed to be given by the Press-Schechter function. The mass density profile 
of dark halos is alternatively taken to be the singular isothermal sphere (SIS), 
the Navarro- Frenk- White (NFW) profile, or the generalized NFW profile. The 
cosmologies being considered include: the Einstein-de Sitter model (SCDM), the 
open model (OCDM), and the flat A- model (LCDM). As expected, we find that 
the lensing probability is extremely sensitive to the mass density profile of lenses 
(dark halos), and somewhat less so to the mean mass density in the universe, and 
the amplitude of primordial fluctuations. NFW halos are very much less effective 
in producing multiple images than SIS halos. For NFW lenses, the SCDM model 
produces lensing events fewer than the OCDM/LCDM models by two orders 
of magnitude. For SIS lenses, the SCDM model produces more lensing events 
with small splitting angles but produces fewer lensing events with large splitting 
angles than the OCDM/LCDM models, which is due to the fact that for large 
mass halos the Press-Schechter function is very sensitive to the amplitude of pri- 
mordial fluctuations. In all cases the difference between the OCDM model and 
the LCDM model is not dramatic. None of these models can completely explain 
the current observations: the SIS models predict too many large splitting lenses, 
while the NFW models predict too few small splitting lenses. Essentially, the 
observed high ratio of small splitting to large splitting lenses is not predicted 
correctly. This indicates that there must be at least two populations of halos in 
the universe: small mass halos with a steep inner density slope and large mass 
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halos with a shallow inner density slope. A combination of SIS and NFW halos 
can reasonably reproduce the current observations, if we choose the mass for the 
transition from SIS to NFW to be ~ 10 13 solar masses as might plausibly occur 
due to baryonic cooling and contraction in lower mass systems. Additionally, 
there is a tendency for CDM models to have too much power on small scales, i.e. 
too much mass concentration. From our sensitivity studies it appears that the 
cures proposed for other apparent difficulties of CDM would help here as well, an 
example being the warm dark matter (WDM) variant which is shown to produce 
large splitting lenses fewer than the corresponding CDM model by one order of 
magnitude. 

Subject headings: cosmology: gravitational lensing — galaxies: clusters: general 
— galaxies: halos 

1. Introduction 

Gravitational lensing directly probes the mass distribution in the universe, so the inves- 
tigation of lensing events of quasars at high redshifts can provide us with important infor- 
mation about cosmology and the matter distribution in the universe. A number of groups 
have attempted to test cosmological models by comparing lensing probabilities predicted 
by various cosmological models with those obtained from observations (Turner, Ostriker, 
& Gott 1984; Narayan & White 1988; Cen et al. 1994; Wambsganss et al. 1995; Maoz et 
al. 1997; Wambsganss, Cen, & Ostriker 1998; Bartelmann et al. 1998; Mortlock & Webster 
2000; Porciani & Madau 2000; Keeton & Madau 2001). In the statistical study of cosmo- 
logical gravitational lensing, two different approaches are used: one is to study the lensing 
probability of splitting angles of multiple images (Turner, Ostriker, & Gott 1984; Narayan 
& White 1988; Cen et al. 1994; Kochanek 1995; Wambsganss et al. 1995; Wambsganss, Cen, 
& Ostriker 1998; Porciani & Madau 2000; Keeton & Madau 2001), the other is to study 
the lensing probability of the length-to-width ratio of arcs formed by gravitational lensing 
(Bartelmann & Weiss 1994; Bartelmann et al. 1998; Cooray 1999; Kaufmann & Straumann 
2000; Meneghetti et al. 2000). The two approaches reflect two different aspects of gravi- 
tational lensing and are complimentary to each other, both of them are directly related to 
observations and deserve detailed investigations. In this paper we focus on the splitting 
angles of multiple images caused by gravitational lensing. The case of length-to-width ratio 
of arcs will be considered in another paper. Results of the work done in testing cosmologi- 
cal scenarios to date are inconclusive with no clear preference found for any of the current 
models, and, overall, one could conclude that none of the studied models provides a good 
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fit to the rapidly improving observations. Perhaps the reason for this indeterminacy is that 
lensing also probes small scale structure. While this is known, its importance has perhaps 
not been appreciated sufficiently Gravitational lensing provides an exquisitely sensitive test 
to the high frequency part of the power spectrum, because it is this part which establishes 
the central profiles of dark matter structures. Since lensing is comparably sensitive to input 
assumptions concerning the properties of lenses and of cosmological models, it is hazardous 
to use observed lensing statistics to draw inferences with regard to cosmology before deter- 
mining the sensitivity to other factors. The primary purpose of this paper is to quantitatively 
assess the sensitivity of lensing expectations to the various input parameters so that this tool 
can be used to greatest effect. 

To determine the probability for gravitational lensing, we need to know the mass density 
profile of lenses (dark halos), the mass function of lenses (which gives the number distribution 
of lenses over masses), the structure and global features of the universe (the universe is flat, 
closed, or open; the matter contents in the universe, etc), and the positions and shapes of the 
source objects (quasars). For studying the splitting angles of multiple images, it is enough to 
assume the source objects are points at high redshifts. In the framework of the Friedmann- 
Robert son- Walker models, there are a number of cosmologies as candidates for describing 
our universe, among which the most popular ones are the Einstein-de Sitter universe, which 
assumes a zero curvature and a zero cosmological constant; the open universe, which assumes 
a negative curvature and a zero cosmological constant; and the flat A-universe, which assumes 
a zero curvature and a positive cosmological constant (Peebles 1993; Ostriker & Steinhardt 
1995). The matter contents in the universe are usually assumed to be dominated by dark 
matter: cold dark matter, warm or hot dark matter, or the mixture of them (Kolb & Turner 
1990; Peebles 1993). Different kinds of dark matter give rise to different mass density profiles 
for dark halos. If the matter contents in the universe are dominated by cold dark matter 
(CDM), the three cosmological models are denoted respectively as SCDM (standard CDM 
= Einstein-de Sitter + CDM), OCDM (open + CDM), and LCDM (cosmological constant 
A + CDM). 

The primordial fluctuations determine the mass distribution and the number density of 
dark halos in the universe. In studying gravitational lensing, the cosmological models, the 
matter contents in them, the power spectrum of the primordial fluctuations, and the source 
object positions are pre-assumed. Then the mass density fluctuations evolve according to 
the Einstein equations, which can be traced with numerical simulations. As time goes on, 
gravitationally bound objects (dark halos) are formed, which play the role of lenses. The 
mass density profile, and the mass function of dark halos in the universe are automatically 
derived from numerical simulations. Then, tracing the trajectories of light rays, we can 
obtain the number of lensing events and thus the lensing probability. This is the basic 
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spirit of numerical simulations for gravitational lensing (Wambsganss et al. 1995). With this 
approach, the mass density profile of lenses and the mass function of lenses appear as the 
results of simulations, they do not need to be pre-assumed. However, numerical simulations 
have their own limitations. First, numerical simulations usually take a great deal of computer 
time, which makes it inconvenient to test several different models at a time. Second, in every 
numerical simulation there are two limits on the spatial scales: the size of box, which is the 
maximum scale and determines the presence or absence of intrinsically rare events, and the 
size of fundamental cell (i.e. the resolution), which is the minimum scale. Correspondingly 
this gives rise the upper and lower bonds on the masses of halos formed in the simulation, 
which causes the result that the number of lensing events with small splitting angles (caused 
by halos with small masses) and large splitting angles (caused by rare halos with large 
masses) are underestimated. This last deficit is seen by comparing the simulation results 
of Wambsganss et al. (1995) with the analytical results of Kochanek (1995), which will also 
be confirmed by the results in this paper. Furthermore, the lensing cross-sections are very 
sensitive to inner most parts of the density profiles which are most subject to problems of 
numerical resolution (Ghigna et al. 2000; Klypin 2000). Ultimately, numerical simulations 
are required to best determine the expected halo properties. The current generation of 
numerical simulations has insufficient dynamical range to simultaneously address accurately 
both large and small scale structures, but forthcoming 1024 3 -particle simulations may remedy 
this defect. For the present, semi- analytical methods may be best adapted to explore the 
sensitivity of the results to the input parameters. 

In the semi-analytical approach, the mass density profile and the mass function of 
lenses are pre-assumed and can be adjusted to fit various dark matter scenarios. Though 
this makes us rely on an assumed density profile and mass function, it has the advantage 
that the calculation of the lensing probability can be handled with semi-analytical methods, 
which takes much less computer time than numerical simulations, so several different models 
can easily be tested at a time. And, in this approach, there are not the limitations in spatial 
scales and mass scales necessitated by numerical simulations; the lensing probability for 
small splitting angles and large splitting angles can be obtained accurately for any assumed 
halo density profile, the most critical feature of which is the steepness of the central density 
profile which may be parameterized by the central slope, —a. The mass density profile is 
often taken to be the singular isothermal sphere (SIS, a — 2) (Gott & Gunn 1974); the 
mass function is usually taken to be the Press-Schechter function (Press & Schechter 1974). 
Though the SIS profile is simple, it does not fit the CDM simulation results well. Navarro, 
Frenk, & White (1996, 1997) have proposed a "universal" density profile (the so-called NFW 
profile, a = 1) which has been shown to fit the simulation results better. A generalized 
NFW (GNFW) profile with somewhat steeper inner slopes intermediate between NFW and 
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SIS appears to provide a still better fit to simulated profiles (Zhao 1996; Subramanian, Cen, 
& Ostriker 2000; Ghigna et al. 2000; Jing & Suto 2000; Wyithe, Turner, & Spergel 2000) for 
which typically a = 1.4 ± 0.1. With N-body simulations of the SCDM scenario (Efstathiou 
1990; Kauffmann & White 1993), or comparison with observations (Bahcall & Cen 1993; 
Girardi et al. 1998; Girardi & Giuricin 2000), the Press-Schechter mass function has been 
shown to describe the mass distribution in the universe quite accurately. 

In this paper we use the semi-analytical approach to investigate the gravitational lensing 
of quasars by foreground dark halos. We work with three kinds of cosmological models: 
SCDM, OCDM, and LCDM models. The mass function of dark halos is taken to be the 
Press-Schechter function. The power spectrum of primordial fluctuations is normalized by the 
cluster abundance constraints (Wang & Steinhardt 1998; Wang et al. 2000). We alternatively 
consider three kinds of density profiles: SIS, NFW, and GNFW profiles, where the last model 
allows us to estimate (by varying the free parameters in the GNFW profile) variants of the 
CDM scenario. For the various models we will calculate the lensing probability of splitting 
angles of multiple images, and test the sensitivity of the lensing probability to the input 
parameters regarding the properties of halos and cosmological models. 

The paper is organized as follows: In section 2 we outline the cosmological models to 
be calculated in the paper, and introduce the power spectrum of primordial fluctuations and 
the Press-Schechter function. In section 3 the most simple case is investigated - SIS halos 
as lenses. We find that, the SCDM model produces fewer lensing events with large splitting 
angles, but more lensing events with small splitting angles than the OCDM model and the 
LCDM model. This is due to the fact that for large splitting angles, which are produced 
by large mass halos, the Press-Schechter function is exponentially sensitive to the amplitude 
of primordial fluctuations and the SCDM model has the smallest amplitude of primordial 
fluctuations; but for small splitting angles, which are produced by small mass halos, the 
Press-Schechter function is not sensitive to the amplitude of primordial fluctuations but the 
mean mass density in the universe plays a more dominant role. In section 4 we consider 
the case for the NFW profile and the GNFW profile. As expected, we find that, the lensing 
probability is very sensitive to the mass density profile of the lenses. The more concentrated 
toward center the density distribution is, the higher lensing probability it produces. Com- 
pared to a SIS halo, a NFW halo is very ineffective in producing multiple images. For the 
NFW profile, the SCDM model produces the least lensing events: the OCDM model and 
the LCDM model produce a comparable number of lensing events which exceed the num- 
ber of lensing events produced by the SCDM model by two orders of magnitudes. This is 
mainly due to two facts: a NFW model in a SCDM universe has the smallest concentration 
parameter, and a SCDM universe has the smallest amplitude of fluctuations. In section 5, 
we compare our predictions with observations. We find that, none of the models can ex- 



- 6- 



plain the current observations: the SIS models predict too many lensing events with large 
splitting angles, while the NFW models predict too few lensing events with small splitting 
angles. Essentially, the observed high ratio of small splitting to large splitting events is 
not predicted correctly. By a combination of them - using the SIS profile for small mass 
halos and the NFW profile for large mass halos, we show that the observation results can 
be reasonably re-produced if we choose the mass at the transition point between the two 
halo populations to be ~ 10 13 solar masses - the mass of halos below which cooling of the 
corresponding baryonic component will lead to concentration of the baryons to the inner 
parts of the mass profile (Rees & Ostriker 1977; Blumenthal et al. 1986; Porciani & Madau 
2000). Our results indicate that the halos in the universe cannot be described by a single 
universal density profile, there exist at least two halo populations: small mass halos with a 
steep inner density slope (the inner parts dominated by the baryonic component), and large 
mass halos with a shallow inner density slope (within which baryons and dark matter have 
essentially the same profiles). In section 6, we draw our conclusions. 

The failure of the SIS, and GNFW (a = 1.5) models at large splitting angles appears 
to be another manifestation of the tendency for CDM models to have too much power on 
small scales, i.e. too much mass concentration. From our sensitivity studies it appears that 
the cures proposed for other difficulties of CDM would help here as well, an example being 
the warm dark matter (WDM) variant, for which a tends to be smaller (Bode, Ostriker, & 
Turok 2000). A model very roughly designed to match a WDM variant has a greatly reduced 
probability of large splitting angles: lower than the corresponding probability for the CDM 
model by one order of magnitude. 



2. Cosmological Models 

The lensing probability is basically determined by three factors: (1) the positions of 
source objects; (2) the positions, the density profile, and the mass function of lenses; (3) the 
cosmological model. In this section we define the cosmological models that will be studied 
in this paper, and introduce the mass function of dark halos. 



2.1. Cosmological functions 



The cosmological models are classified with parameters Q m , and Qr, which are 
respectively defined by (Peebles 1993) 

= ^ , = — , = 7— ^ 2 ' (1) 

Pcrit.o o7rLrp cr it,o o7rLrp cr it,o0 
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where po is today's average mass density (dark matter + baryonic matter + radiation) in the 
universe, G is the gravitational constant, A is the cosmological constant, a is today's scale 
of the universe, k = 0, ±1 is the spatial curvature of the universe, p C rit,o is today's critical 
mass density of the universe which is defined by 

Pcrit,o = |§ « 1-88 x KT 29 h 2 g cm" 3 , (2) 

where Hq = lOO/i km s _1 Mpc^ 1 is the Hubble constant. The Einstein equations tell us 
that Q m + £7a + flu = 1. Thus, among Q m , Q^, and only two of them are independent 
parameters. In this paper we assume that the matter in the universe is dominated by cold 
dark matter (CDM), and we consider three kinds of universes: 

1. Einstein-de Sitter model (SCDM). In this model, Q m — 1, fl\ — = 0, and k = 0. 
The linear growth function, which describes the growth speed of the linear perturbation in 
the universe (Peebles 1980), is 

fW~. (3) 

where z is the cosmic redshift. The linear growth function is normalized to be unity at 
present epoch: V(z = 0) = 1. 

2. Open model (OCDM). In an open universe, < Vl m < 1, £l\ — 0, < — l — VL m < 
1, and k = — 1. The linear growth function is 

where 

Fl ( X ) = 1 + I + 3 -^±f^ In [(1 + xY<* - x^] . (5) 

Fi is fitted by 

1.5 ^ a6 



with an error < 1.3% for 10~ 5 < x < 10. 

3. A-model (LCDM). A A-universe has < Q, m < 1, < fi A = 1 - Q, m < 1, Q, R = 0, 
and k = 0. The linear growth function is 



1/3 



(2n A /n m ) 

1 + 2 
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where 



Fo(x) = 



X 



1/2 



u 



u 3 + 2 



3/2 



F 2 is fitted by 



F 2 {x) « 0.358 [l - (1 + 0.23a; 2 - 2 ) 1-26 
with an error < 2% for < x < 10 4 . 



du . 



1/2.2 



(8) 



(9) 



For all the three cosmological models, the proper cosmological distance from an object 
at redshift z\ to an object at redshift z 2 along the same line of sight is (Peebles 1993) 



D (z 1 ,z 2 ) 



Z2 



(l + z)- 1 dz 



H J Z1 y/n m (i + zf + n R (i + z y + n A ' 



(10) 



where c is the speed of light. In studying gravitational lensing, another useful distance is 
the angular-diameter distance. The angular-diameter distance from an object at redshift z\ 
to an object at redshift z 2 along the same line of sight is (Bartelmann & Schneider 2001) 



-1/2 



D A ( Zl ,z 2 ) 



c/H 
l + z 2 



sin 



w 



f2 R 1//2 sinh 



R\ 



1 1/2 



W 



where 



w (zi,z 2 ) = 



n 1 / 2 
\l R ' w 



dz 



(*=1), 
(k = 0) , 
(k = -l), 



(11) 



21 



y/n m {i + z)* + n R (i + z) 2 + n A 



(12) 



Note that D(z u z 2 ) = D (z 2 , z±) and D(z u z 2 ) = D(0,z 2 ) -D(0,zi), but D A (^i,2 2 ) 7^ 
L> A (22, zi) and D A (z u z 2 ) ^ L> A (0, z 2 ) - -D A (0, z x ). But we have 



(1 + z 2 ) D A ( Zl , z 2 ) = (1 + ax) £ A (z 2 , Zl ) 



(13) 



2.2. Mass fluctuations 

Dark halos of galaxies and galaxy clusters are formed from primordial fluctuations of 
matter in the early universe. For a fluctuation with a power spectrum P k = \5k\ 2 , its variance 
on a scale of comoving radius r is (Kolb & Turner 1990) 

— \ =W 3 J q ^k 2 dkP k W 2 (kr), (14) 
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where W(kr) is the Fourier transformation of a window function. For a top-hat window 
function 



W(kr) 

for a Gaussian window function 



sin kr cos kr 



(kr) 3 (kr)' 



W(kr) = exp[-k 2 r 2 /2] 



(15) 



(16) 



In this paper we use the top-hat window function. We compute the CDM power spectrum 
using the fitting formulae given by Eisenstein & Hu (1999) 



— P(k,z) = 6' H [ — j T\k)V(z)\ 



3+n 



(17) 



where is the amplitude of perturbations on the horizon scale today, n is the initial power 
spectrum index, V(z) is the liner growth function, and 

(18) 



L + C£„ ' 



where 



ln(e + 1.84g e// ) , q eff 



k 



VL m h? Mpc 



-l ' 



C = 14.4 + 



325 



1 + 60.5^// ' 

We normalize the power spectrum to er| = A 2 (z = 0, r = 8/i _1 Mpc), then 



(19) 
(20) 



>H 



/ \ 3+n 

JT f (it) T2 (*) ^ 2 M 



1/2 



(21) 



The value of a 8 can be estimated from the cluster abundance constraints (Wang & Steinhardt 
1998; Wang et al. 2000) 



a 8 f^ « 0.5 , 



(22) 



where 7 « 0.43 + 0.33fi m for LCDM and 7 « 0.33 + 0.35fi m for OCDM. For a scale- 
invariant spectrum predicted by the inflation theory and being consistent with observed 
cosmic microwave background (CMB) fluctuations, we have n — 1 (Kolb & Turner 1990). 
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2.3. Mass function of dark halos 

Assume the primordial density fluctuations are Gaussian. Then, according to the Press- 
Schechter theory (Press & Schechter 1974), the comoving number density of dark halos 
formed by redshift z with mass in the range (M, M + dM) is 

n(M,z)dM = ^f(M,z)dM, (23) 

where p = ^mPcrit.o is the present mean mass density in the universe, f(M, z) is the Press- 
Schechter function 



[2 SJz) din A 



2A 2 



(24) 



where A 2 = A 2 (M, z = 0) is the present variance of the fluctuations in a sphere containing 
a mean mass M, and S c (z) is the density threshold for spherical collapse by redshift z. The 
redshift- dependent density threshold has been calculated by many people (Lacey & Cole 
1993; Kochanek 1995; Eke, Cole, & Frenk 1996, and references therein), here we use the 
approximation (Navarro, Frenk, & White 1997) 

where 

o/ x = _P_ _ Am(l+*) 3 (26] 

[, - Pciit n m (i + z) 3 + n R (i + z y + n A > 1 ) 

T>(z) is the linear growth function normalized to unity at z — 0, and 

{1, if Q m = 1 and A = 

fi - 0185 , if tt m < 1 and A = . (27) 
fi - 0055 , if Vt m + n A = 1 

To a very good approximation 5 c (z)V(z) ~ 1.69 for all cosmological models with Q m > 0.1. 
In equation (24), r is related to M by 

M = yftr 3 = — n m Pcrit.o r 3 , (28) 



i.e. 



1 M \ 1/3 
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Though it is obtained from very simple considerations, the Press-Schechter function has 
been shown to be in remarkable agreement with N-body simulations for the standard cold 
dark matter scenario (Efstathiou 1990; Kauffmann & White 1993) and observations (Bahcall 
& Cen 1993; Girardi et al. 1998; Girardi & Giuricin 2000). 

The fluctuation variance A 2 decreases with increasing mass M . Thus, from equation 
(24), f(M, z) decreases exponentially with increasing M. Because of this, a crude knowledge 
of f(M, z) at the large mass end gives a strict constraint on cosmological parameters (Chiu, 
Ostriker, & Strauss 1998). To see this, let us write the Press-Schechter function in the form 



f(M, z, a 8 ) = : exp 



A(M, z) 



(30) 



where we have used the fact that A oc <r 8 . From equation (30), we have 



f 



^8 



2 



(31) 



So, if Aja\ >> 1/2 - which is true for large M - a large error in / corresponds to a small 
error in a 8 . As a simple example, consider a power law spectrum P k oc k ni , then we have 



A 2 (M) = al 



1.189 



1 



M 



1/3' 



-ni-3 



(32) 



and 



2A 



. "1+3 



oc M— D(z)- 2 a 8 



(33) 



For large M and large z, a large 5f / f corresponds to a small Sag/ag if rii > —3. 

In the upper panel of Fig. 1 we plot the Press-Schechter function / against the velocity 
dispersion a v of dark halos at z = in a SCDM cosmology, where the CDM power spectrum 
is given by equation (17). The velocity dispersion of a dark halo is defined by 



/ GM\ 
V2r 2 oo/ 



1/2 



(34) 



where r2oo is the radius of a sphere around a dark halo within which the average mass density 
is 200 times the critical mean mass density of the universe, i.e. 



M 



4?r„3 
3 '200 



200p, 



crit i 



(35) 
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and M = M 2 oo is the mass within that sphere 



'"200 



M = M 200 = 4vr / pr 2 dr . (36) 



o 



From equation (34) and equation (35), we have 
.3 / q \ 1/2 



M 



G \100nGp crit/ 

0.656 x lO 15 / i ~ 1 M crj 1000 [fi m (l + zf + + zf + ft A ] ~ 1/2 , (37) 



where M Q is the solar mass, cr^iooo = ioook m /s ' an< ^ we have used 

Pcrit(^) = Pcrit,o (1 + zf = p criti0 [Vl m (l + zf + n R (l + zf + n A ] . (38) 

In the upper panel of Fig. 1, the solid curve is for the case of a 8 = 0.5, the dashed curve 
is for the case of a 8 = 0.6. It clearly shows that with the same amount of change in a 8 , a 
larger change in / happens at the end of large a v . In the lower panel of Fig. 1, we plot 
£ = /JJ^ against a v for a SCDM cosmology with z = and a 8 = 0.5. Again, it shows 
that for large a v a poor knowledge in / gives a good estimation of a%. This is true because 
the number density of halos with high velocity dispersion (which is true for rich clusters of 
galaxies) is very sensitive to the value of <r 8 . Thus, since lensing is produced primarily by the 
highest velocity dispersion clusters, the observed number of lenses will sensitively constrain 
the value of cr 8 . 



3. Lensing by a Singular Isothermal Sphere 

3.1. Singular isothermal sphere as a lens 

The most simple model for a lens is a singular isothermal sphere (SIS) with a mass 
density 

M = £ah< (39) 

where a v is the velocity dispersion (Gott & Gunn 1974; Turner, Ostriker, & Gott 1984; 
Schneider, Ehlers, & Falco 1992). Though it is simple, this model can describe the flat 
rotation curves of galaxies and many basic features of gravitational lensing. Though the 
gravitational lensing by a SIS has been well studied in the literature (Gott & Gunn 1974; 
Turner, Ostriker, & Gott 1984; Narayan & White 1988; Kochanek 1995), we present it here 
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since we are using an updated CDM power spectrum and we want to compare its results 
with those for the NFW and the GNFW cases. 

The surface mass density of the SIS is 

e «>=iB' (4o) 

— # — # 

where £ = £ is the position vector in the lens plane. Choose the length scales in the lens 
plane and the source plane to be respectively 



where Dg is the angular-diameter distance from the observer to the source object, is the 
angular-diameter distance from the observer to the lens object, D^ s is the angular-diameter 
distance from the lens to the source object. Remember that D^ s ^ Dg — D^. Then the 
position vector of a point in the lens plane can be written as £ = x£o, the position vector of 
a point in the source plane can be written as if = yi] . The lensing equation is 

mix) 

y = x- ^-^ , (42) 

x 



where 



m(x) = 2 [ X ^^x'dx' , (43) 



o 



where the critical surface mass density S cr , defined by (Turner, Ostriker, & Gott 1984) 



c 2 



Scr " 4ttG D£D£ S ' (44) 

is of the order of the surface density of the universe. 

Inserting equation (40) and equation (44) into equation (43), then into equation (42), 
we obtain S/S cr = (2|a:|)~ 1 and 

y = x- — . (45) 

x 

Double images are formed if and only if \y\ < 1, i.e. \x\ < 1. The separation between the 
two images is 



Ax = 2, 



(46) 
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thus the splitting angle is 

AS = J^Ax = 8, 5& = 0.9613' < 1000 5& . (47) 

5 5 

Using equation (37), we have 

_ 8vr /lOOvr V/ 3 _ / M \ 2 / 3 / p crit V/ 3 

A " " ^ l"r G M Pcrit J -of ~ 1271 ~Df Wsh-iMj l^J • (48) 

The cross-section (defined in the lens plane) for forming two images with splitting angle 
Ad > A9 is 

a = tt£ 2 (A6 - A9 ) = tt£ 2 i? (M — M ) , (49) 
where i? is the step function 



and 



/ 3 y/ 2 C 3 (Ae^^fDlV' 2 



3/2 , nA x 3/2 , \ -1/2 



0.6975 x 10^ (*») (^) (-^j . ,5!) 



3.2. Lensing probability by a singular isothermal sphere 

The integral lensing probability for a source object at redshift z s is (Schneider, Ehlers, 
& Falco 1992) 

fcfe, (52) 



o 



where z is the redshift of lenses, and 
dP dD L 



oc 



fi(M, z)a(M, z)dM 



dz dz Jq 

dD f°° 1 

= Pcrit,o^ m (l + ^) 3 ^y o —f(M,z)a(M,z)dM, (53) 



- 15 - 



where n(M, z)dM = n(M,z){\ + z) 3 dM is the physical number density of dark halos of 
masses between M and M + dM, n(M, z)dM is the comoving number density of dark halos 
of masses between M and M + dM, a(M, z) is the lensing cross-section for a dark halo of 
mass M at redshift z, f(M,z) is the mass function of dark halos. Here Dl is the proper 
distance from the observer to the lens object. 

If we use 1O 15 /i _1 M as the unit of M, then the unit of f(M, z) is (1O 15 /i -1 M ) -1 . Define 
Mi 5 = M/(lO 15 /i _1 M ). Let us use c/if = 2997.9 /i _1 Mpc as the unit for cosmological 
distances, and ft, _1 Mpc as the unit for local lengths (thus the unit of a is /i~ 2 Mpc 2 ), then 



d 4- = 0.8321 n m (l + z) 3 j- , , 

u2 \c/H J J Q 



f(M 15 ,z) 



v(M 15 ,z) 
lh~ 2 Mpc 2 



dM 



15 



Mi 



15 



(54) 



Inserting equation (49) and equation (41) into equation (53), we have for the SIS case 



^P(> A6 ) = 167T 3 P cnt,0fi m (l + Zf 



dD± (D£D£ S 

dz 



u s 



M 



±f(M,z)(^ydM, (55) 



where a v /c is related to M by the reverse of equation (37), i.e 

1/3 



a v 
c 



[n m {i + zf + n R {i + zf + n K ] 1/& 

3.836 x 10 _3 M 1 1 5 /3 [n m (l + zf + Q R (1 + zf + Q A ] 1/6 



where 



_ / 3c 6 N 1/2 
Ml = Vl007rG 3 p criti0 , 



1.772 x 10 22 /i _1 M m . 



o ■ 



Or, from equations (49), (41), and (54), we have 



|P(>A W = 0.4593 Ml + -fL(*L 



X 



/■OO 

/ /(M 15 , 

J M 



z) 



( ^ 



V10 3 km/sy Mis ' 



dM 



15 



where = D^D£ s /D§ . 



(56) 



(57) 



(58) 



With equations (52), (56), and (58), we can calculate dP(> A9 )/dz and P(> A9 ) 
for a given cosmology (specified by Q m , and the Hubble constant h), a given primary 
perturbation (specified by a 8 ), and a given position of the source object (z s ), assuming the 
CDM power spectrum is given by equation (17) and the mass function of dark halos is given 
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by equation (24). In our numerical calculations we assume h = 0.7, z s = 1.5, and a 8 is 
related to Q m by equation (22). In Fig. 2 and Fig. 3 we show the results for three typical 
cosmological models: SCDM with Q m = 1 and a 8 = 0.5; OCDM with Q m = 0.3, Qa — 0, 
and a 8 = 0.85; and LCDM with Q m = 0.3, Q A = 0.7, and a 8 = 0.95. 

Fig. 2 shows the differential lensing probability dP/dz for A9 > A9 = 5". The maxi- 
mum of dP/dz is at z « 0.3 for SCDM, at z « 0.4 for OCDM and LCDM. Median expected 
redshifts for the three models are respectively 0.349, 0.455, and 0.461. This is consistent 
with the demonstration that the differential lensing probability peaks at intermediate red- 
shifts (i.e. roughly the half way to the source object; Turner, Ostriker, & Gott 1984) and 
with the scenario that structures form later in the SCDM model than in OCDM and LCDM 
models (Longair 1998; Bartelmann et al. 1998). As is well known structure forms earlier in 
low Q m universes, so the relative probability of a lens being at z — 0.1 rather than 1.0 is 
orders of magnitude higher in SCDM than it is in OCDM/LCDM. Fig. 3 shows the integral 
probability P(> A9) as a function of the splitting angle A9. For small A9, the SCDM model 
produces more lensing events than the OCDM and the LCDM models, which is a manifes- 
tation of the fact that the SCDM universe has a higher mean mass density than the OCDM 
and the LCDM universes. For large A9, the SCDM model produces less lensing events than 
the OCDM and the LCDM models, which is caused by the fact that at large mass end the 
Press-Schechter function is very sensitive to a 8 (see Fig. 1) while the SCDM model has 
a much smaller a 8 than the OCDM and the LCDM models. The difference in the lensing 
probabilities for the OCDM model and the LCDM model is mainly due to the effect of the 
cosmological constant: the cosmological constant tends to increase the lensing probability 
(Turner 1990), primarily because there is a greater metric distance and thus more potential 
lenses between an observer and a source with a given redshift in A cosmologies. 



4. Lensing by a Generalized NFW Model 

4.1. The generalized NFW model as a lens 

The SIS model is simple and can describe many observational features of galaxies and 
galaxy clusters. However, the SIS does not fit well the density profiles predicted by N-body 
numerical simulations (Frenk et al. 1985, 1988; Quinn, Salmon, & Zurek 1986; Efstathiou et 
al. 1988; Zurek, Quinn, & Salmon 1988; Warren et al. 1992; Crone, Evrard, & Richstone 
1994). Navarro, Frenk, & White (1996, 1997) proposed a "universal" mass density profile 
for dark halos: the so-called NFW profile, which fits the simulation results better than the 
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SIS. The NFW profile is 



3 

PNFw(r) = . Ps _[ s V2 , (59) 
r[r + r s ) 2 



where p s and r s are constants. At large radii (r ^> r s ) Pnfw ~ r_3 > & t small radii (r <C 
r s) Pnfw ~ both are different from the SIS. At intermediate radii the NFW profile 
resembles the SIS. Though Navarro, Frenk, & White demonstrated that over two decades in 
radius their profile accurately fits halos with mass spanning about four orders of magnitude 
ranging from dwarf galaxy halos to those of rich galaxy clusters, with higher resolution 
simulations Jing & Suto (2000) have argued that the NFW profile is not correct on small 
scales. Furthermore, Subramanian, Cen, & Ostriker (2000) and Ghigna et al. (2000) have 
emphasized the dependence of the inner slope of halo mass density on the form of the power 
spectrum of the primordial fluctuations. Therefore, here we consider a generalized NFW 
(GNFW) profile (Zhao 1996) 

where a (0 < a < 3) is a new constant parameter. The NFW profile is a specific case of the 
GNFW profile when a — 1. Obviously, at large radii (r 3> r s ) the GNFW profile has the 
same behavior as the NFW profile. But on small scales (r r s ) they are different unless 
a — 1. When a — 2, the GNFW profile resembles an SIS on small scales but resembles the 
NFW profile on large scales. 

The surface mass density for the GNFW profile is 

E(x) = 2p s r s {x 2 + z 2 )~ a/2 \{x 2 + z 2 ) 1/2 + l] dz , (61) 

where x = \x\ and x = £/r s , £ is the position vector in the lens plane. Inserting equation 
(61) into equation (42) and equation (43), we obtain the lensing equation for a GNFW halo 

y = x-p s 9 -^l, (62) 

x 

where y — \y\, ff — y r s Dg / is the position vector in the source plane, and 

g(x)= udu (u 2 + z 2 y a/2 \(u 2 + z 2 ) 1/2 + l\ dz, (63) 
Jo Jo L 



and 



Ps = , (64) 



J cr 



-18- 



where S cr is the critical surface mass density defined by equation (44). Note, here we use 
£o = r s as the length unit in the lens plane, r] = r s Dg /D^ as the length unit in the source 
plane. The dimensionless parameter fi s summarizes the ability for a GNFW halo to produce 
multiple images. Multiple images are formed if and only if \y\ < y cr , where y CT = — y(x CT ), 
x cr > is determined by dy/dx = (see Fig. 4). Let us consider the splitting angle A9 
between the two outside images when more than two images are formed. For \y\ < y cr , there 
are three real roots of equation (62): x\ > x 2 > x 3 , then AO oc Ax = x 1 — x 3 . In general, the 
value of X\ — x 3 (and thus the value of A9) is insensitive to the value of y when \y\ < y CT 1 . 
So, we have 

Ax (y) Ax (y = 0) = 2x , for |y| < y cr , (65) 

where x is the positive root of y(x) = 0. Then, for a GNFW lens at redshift z, the cross- 
section in the lens plane for forming multiple images with A# > A9 is 

<j (> A#0i M, z) ~ ny 2 CI r 2 s (AO - A6 ) . (66) 

The splitting angle A9 is given by 

M = ^ AX ^ 2J P- (67) 



For a — 1 (the NFW case) and a = 2 (the modified SIS case), g(x) defined by equation 
(63) can be worked out analytically. For a — 1, we have 



x 



== arctan \/ x 2 — 1 (x > 1) 



g(x) = In - + < -y=f arctanhVl - x 2 (0 < x < 1) , (68) 

! 1 (.r 1 ) . 



For a = 2, we have 



-y/x 2 — 1 arctan \/x 2 — 1 (x > 1) , 



g(x) = In ^ + ^-x + <J VI - x 2 arctanhVl - x 2 (0 < x < 1) , (69) 



2 2 



1 (x = 1) . 



Equation (68) has also been obtained by Bartelmann (1996). For other values of a, g(x) has 
to be worked out numerically. 



1 See Schneider, Ehlers, & Falco (1992). For the extreme case of a singular isothermal sphere, Ax is 
exactly independent the position of the source object in the source plane when double images are formed, 
see equation [46]. 
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4.2. Determination of p s and r s 

For a halo with a GNFW profile, its mass diverges logarithmically as r — > oo. So as 
usual, we define the mass of a halo to be the mass within r 2 oo 

rr200 

M = 4tt pr 2 dr = 4irp s r 3 J( Cl ), (70) 
Jo 

where c\ = r 2 oo/r s is the concentration parameter and 

For 1 < a < 2, /(ci) can be worked out analytically 

f Ind + d)-^ ( tt = l), 

/(d) =^ Ml + ci) (a = 2), (72) 

[ |-I 2 Fi(3-a,3-a;4-a;-ci) (1< a < 2) , 

where 2 Fi is the hypergeometric function. From equation (35) and equation (70), we obtain 

200 c 3 200 c 3 

Ps = Pcrit — = Pcrit.O [fi m (1 + ^) 3 + + zf + fi A ] — ~77~~T , (73) 

3 /(Ci) 3 /(Ci) 



1 / 3M \ 1/3 1.626 Mil 3 l-u* 

= — I TTTTT I = 77TT k MpC . (74) 



and 



Ci V8007T Pcrit ; " Ci [(7 m (l + 2 )3 + (7 i? (l + 2 )2 + (7 A ] 1 / 3 

If we know the value of the concentration parameter ci, with equation (73) and (74) we 
can determine p s and r s for any halo of mass M in any cosmology. Interestingly, p s does 
not depend on the mass of the halo. To determine the value of C\ is not easy, and different 
methods give different results (Navarro, Frenk, & White 1997; Bartelmann et al. 1998). 
Various simulations suggest that c\(z) oc (1 + z)^ 1 (Navarro & Steinmetz 2000; Bullock et 
al. 2001). Thus, here we assume that C\{z) — C\(z — 0)/(l + z), and we try to infer c±(z = 0) 
from simulation results. We take the values of C\(z = 0) for the NFW profile (i.e. the 
a = 1 case) from Bartelmann et al. (1998)'s simulation results, then obtain the values of 
c\(z = 0) for other cases (i.e. a > 1) by referencing the values for the NFW case. To do 
so, let us assume that fitting a dark halo with different density profiles gives the same ratio 
r] = ri/ 2 /r 2 oo, where ry 2 is defined by M(r < ri/ 2 ) = \M(r < r 20 o)- Then we obtain a 
relation 

T ' ci(z=0) x 2 dx _ i r ci(z=0) x 2 dx 

x a (l + x) 3 ~ a ~ 2 J Q x a (l + xf~ a ' ^ ' 
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In equation (75), i] is the same for all values of a. So, with the known value of c\ for the 
a = 1 case - let us denote it with c , we can solve 77 = r)(c ) by equation (75). Then, since 
77 is assumed to be the same for all the values of a, we can solve c\(z = 0) = c\(a, z = 0) by 
equation (75) for any 1 < a < 2 with the rj just solved. Then, we calculate C\ at any redshift 
with a(z) = ci(z = 0)/(l + z). 



According to Bartelmann et al. (1998), for the NFW profile, we choose 



ci(z = 0) 



5, forSCDM 
9, for OCDM 
7, for LCDM 



(76) 



Then, according to the procedure described above, we obtain 



ci(z = 0) 



for a GNFW profile with a — 1.5, and 



ci(z = 0) 



2.7, for SCDM 
5.3, for OCDM 
4.0, for LCDM 



0.58, for SCDM 
1.8, for OCDM 
1.2, for LCDM 



(77) 



(78) 



for a GNFW profile with a = 2. 

A critical parameter determining gravitational lensing is the surface mass density. A 
mass concentration with a central surface density larger than the surface density of the 
universe, E cr oc cH /G (see eq. [44]), can produce multiple images (Turner, Ostriker, & 
Gott 1984). While it is true that for a > 1 the surface density is divergent as r — > 
and thus, formally, all halos can produce multiple images, there may be very little mass 
contained within the E = E cr contour. So, large splittings will not be common unless the 
surface density at the half- mass point is near E cr . The mean surface density at the half- mass 
radius is 



J l/2 



M 
2^ 



8007rp c 



2/3 



1 



3M / ?r 

0.0126 ^ M\i z h [Q m (l + zf + n R (i + zf + Q A ] S-(ci) , 



(79) 



where S(ci) = rj(ci)~ 2 . For a fixed cosmology and a fixed mass Ex/2 is proportional to S(ci). 
We plot this function in Fig. 5. The steep dependence of S(ci) (oc S/E cr ) on c\ indicates 
that, for any assumed a, more concentrated halos are much more effective in lensing. 
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Several authors have noted that the large concentrations indicated by high resolution 
N-body simulations of the CDM scenario (Navarro, Frenk, & White 1996, 1997; Moore et al. 
1998) may be inconsistent with a variety of observations including the inner rotation curves of 
galaxies (Moore 1994; Flores & Primack 1994; De Block & McGaugh 1997; Tyson, Kochanski, 
& Dell'antonio 1998; Spergel & Steinhardt 2000). Here we note that if we alter, for whatever 
reasons, the concentration parameter from ten to five we will lower the characteristic half- 
mass surface density by a factor of ~ 1.5 reducing greatly the fraction of mass in the universe 
contained in halos with a mean surface density exceeding the critical surface density, and thus 
reducing the probability of lensing by a large factor. This will be confirmed by the results 
in the next subsection. Recently, Bode, Ostriker, & Turok (2000) have performed detailed 
high-resolution N-body simulations of the warm dark matter (WDM) scenario to determine 
if this variant can successfully address the putative difficulties of the CDM paradigm. They 
find a significant decrease in concentration in the WDM scenario and also a decrease in the 
best fit value of a. 



where equations (44) and (64) have been used. With equations (62), (67), (73), (74), and 
(80), we can solve A9 = A9(M,z) for any 1 < a < 2. In Fig. 6 we show A9 = A9(M) 
produced by a NFW lens at z — 0.3, where the concentration parameters are given by 
equation (76). The source object is assumed to be at z s = 1.5. For comparison, we also 
show the splitting angle produced by a SIS lens at the same redshift (thin lines in Fig. 6). 
The difference in the splitting angles produced by the two different mass density profiles 
is dramatic for small mass lenses: the splitting produced by a NFW lens is greatly shifted 
towards small angles. To produce the same small splitting angle a NFW halo requires more 
mass than a SIS halo. And, for the NFW case the difference in the results for different 
cosmological models are important - especially for low mass lenses, but for the SIS case the 
difference in the results for different cosmological models are unimportant. 

The splitting angle A9 increases monotonically with increasing M, so $(A9 — A9 ) = 
-&(M - M ), where M is obtained by solving A9(M ,z) = A9 . Then, from equation (54) 
and equation (66), the differential lensing probability by a GNFW profile is 



4.3. Lensing probability by a generalized NFW model 



Once p s and r s are determined, we can calculate ji s with 




(80) 



d 

dz 
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/ f(M 15 ,z)y 2 cl 



r_s \ 2 dM 15 



(81) 



J M 



lh- l Mpc) Mi 5 ' 



The integral lensing probability is calculated with equation (52) and equation (81). 

We have calculated the lensing probabilities for different cosmologies and different 
GNFW halos. Our results are summarized in Fig. 7 - Fig. 10, where the source object 
is assumed to be at z s = 1.5, the Hubble constant is taken to be h = 0.7, cr$ and Q m are 
constrained by equation (22). 

Fig. 7 shows the differential lensing probability dP/dz for A9 > 5" as a function of 
lenses' redshift z. The three cosmological models are: SCDM with fl m = 1 and a 8 = 0.5; 
OCDM with tt m = 0.3, tt A = 0, and a 8 = 0.85; LCDM with Q m = 0.3, Q A = 0.7, and 
<r 8 = 0.95. Three cases with different values of a are shown: a — 1 (i.e. the NFW case), 1.5, 
and 2. The concentration parameters of halos are determined with the procedure described 
in the last subsection, i.e. they are given by equation (76), equation (77), and equation 
(78). Similar to the case for SIS halos (see Fig. 2) and independent of the parameter a of 
GNFW halos, the differential lensing probability peaks at intermediate redshifts: at z ~ 0.3 
for SCDM, at z ~ 0.4 for OCDM and LCDM. Fig. 8 shows the integral lensing probability 
P{> Ad) as a function of the splitting angle A6, for the same models in Fig. 7. As a 
decreases (i.e. as the halo central density becomes more shallow), Fig. 8 shows, the lensing 
probability drops quickly. NFW halos are least efficient in producing multiple images among 
halos with 1 < a < 2. Furthermore, as a decreases, the difference between the SCDM 
cosmological model and the OCDM/LCDM model becomes more prominent: for the NFW 
case the lensing probability for the SCDM model is lower than the lensing probability for 
the OCDM/LCDM model by two orders of magnitude, though the difference between the 
OCDM model and the LCDM model is not so big. This is due to the fact that the lensing 
probability for the NFW case is extremely sensitive to the concentration parameter of halos 
and the mass fluctuation <7 8 , a halo in a SCDM universe has the smallest concentration 
parameter c\ (see eq. [76]), and the SCDM model has the smallest <t 8 . Fig. 8 shows one 
of our most important results. If a is the same for all halos, then one expects that, in the 
OCDM/LCDM case, the number of 10" splittings observed to be only slightly less than the 
number of 1" splittings observed. 

The sensitivity of the lensing probability to the concentration parameter is shown in 
Fig. 9, where the models are the same in Fig. 7 and Fig. 8 except that in Fig. 9 we allow 
the concentration parameter to vary from 2 to 15. Fig. 9 shows that the sensitivity of the 
lensing probability to the concentration parameter increases quickly as a and C\ decrease. 
Fig. 10 shows the dependence of the lensing probability P(> 5") on the cosmic mass density 
Q m when Q m and erg are correlated through equation (22), the observed cluster number 
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constraints. Though the correlation in equation (22) is confirmed by the statistics of cluster 
abundances, we see some breakup of this relation in lensing statistics. The breakup is most 
prominent for the NFW case. 



5. Comparison with Observations 

To compare with observations we must consider the effect of magnification bias (Turner, 
Ostriker, & Gott 1984; Kochanek 1995; Bartelmann & Schneider 2001). According to Turner, 
Ostriker, & Gott (1984), the factor B by which lensed objects at redshift z s will be overrep- 
resented in any particular observed sample may be written as 

_ J-S(f)JZA^P(A)N Zs (f/A)dAdf 

JT s(f)N Zs (f)df ' [ * Z) 

where / is the observable flux of source objects, S(f) is the selection function, A m is the 
minimum total flux amplification for a multiple imaged source object, P(A) is the probability 
density for a greater amplification A, and N Za (f)df is the number of source objects in the 
sky at at redshift z s with unlensed flux lying between / and / + df. If the sample of source 
objects has a flux distribution with a single power-law N Za oc / _/3 , the calculation of the 
magnification bias factor becomes extremely simple 

roc 

B= / A f3 ~ 1 P(A)dA , (83) 



which is independent of the selection function. Quite generally P(A) is given by P(A) = 
1A 2 m A- z for A > A m (Turner, Ostriker, & Gott 1984; Schneider, Ehlers, & Falco 1992). 
Then, if the sample of source objects has a power-law flux distribution, the magnification 
bias is simply 

B-^Ai-K (84) 

For a SIS lens the total magnification is A = 2/|y| so A m = 2 since multiple images 
are formed only if y < 1 (Turner, Ostriker, & Gott 1984; Schneider, Ehlers, & Falco 1992). 
For a NFW or GNFW lens, the calculation of the total magnification is straightforward but 
complex, and to determine its minimum is not easy. Here we estimate the minimum of the 
total magnification for a NFW/GNFW lens with 

A m w ^ , (85) 
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where xq is the positive root of y(x) = 0, y cr = —y(x CT ) and x cr is the positive root of 
dy/dx = (see subsection 4.1 and Fig. 4). Note equation (85) is accurate for a SIS lens. 
Using equation (66) and equation (67), equation (85) can be written as 



where a is the lensing cross-section. 

For a SIS lens A m = 2 which is a constant number, thus B is independent of the redshift 
and mass of the lens and the observable lensing probability P G bs is related to the intrinsic 
lensing probability P simply by P obs = BP (Turner, Ostriker, & Gott 1984). But for a 
NFW/GNFW lens, A m and B depend on both the redshift and the mass of the lens, thus 
Pobs is related to P by an integration 



where the intrinsic lensing probability P is given by equation (52) and equation (53). The 
average magnification bias B = P bs/P is then a function of A0. 

Currently the largest uniformly selected sample of gravitational lensing system is the 
Cosmic Lens All-Sky Survey (CLASS; Browne & Myers 2000). The sample comprises 11685 
flat-spectrum radio sources whose flux distribution is given by a power-law N Zs (f) oc J" 2 ' 1 
(Helbig 2000; Rusin & Tegmark 2000). The redshift distribution of the CLASS sample is 
not known, but Marlow et al. (2000) has reported redshifts for a small subsample of 42 
sources. They have found a mean source redshift of (z s ) = 1.27, which is comparable to 
that found in other radio surveys at comparable fluxes. To date a total of 18 multiple image 
gravitational lenses have been discovered in the combined JVAS 2 and CLASS sample, all 
have image separations A9 < 3" (Browne & Meyer 2000; Helbig 2000). An explicit search 
for lenses with images separations 6" < A9 < 15" has found no lenses (Phillips et al. 2000). 
Among the 18 discovered lens systems, one (B2114+022) is questionable (Helbig 2000) so 
we exclude it from our analyses. The remaining 17 lens systems are shown in Fig. 11 as a 
histogram. 

To compare our results with the CLASS survey, we assume all source objects in the 
sample are at z s = 1.27. For the CLASS sample f3 ~ 2.1 (Rusin & Tegmark 2000), so from 
equation (84) we have 




(86) 




(87) 



B « 2.22 A, 



1.1 



(88) 



rn 



2 Jodrell-VLA Astrometric Survey, see King et al. (1999) and references therein. 
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For SIS lenses A m = 2, then we have B « 4.76. For NFW/GNFW lenses we calculate A m 
and B using equation (86) and equation (88). Then we calculate the observable lensing 
probability using equation (87). The results for SIS lenses and NFW lenses are shown in 
Fig. 11 alternatively for SCDM with Vt m = 1 and a 8 = 0.5; OCDM with Vt m = 0.3, fi A = 0, 
and as = 0.85; LCDM with Q m = 0.3, Q\ = 0.7, and og — 0.95. Again, we assume the 
Hubble constant h = 0.7. The concentration parameters for the NFW lenses are c\ — 5 
for SCDM, 9 for OCDM, and 7 for LCDM. For comparison, we have also shown the result 
for a LWDM model with a dashed-dotted curve, which has the same parameters as the 
LCDM model except that for the LWDM we assume a power-law spectrum P k oc k~ 2 and 
take a smaller concentration parameter c\ = 3.5. The null results for the JVAS/CLASS 
survey for 6" < A6 < 15" is shown with a horizontal straight line with a downward arrow 
indicating that is an upper limit. Fig. 11 clearly shows that, all SIS models predict too many 
lenses for large splitting angles. While for large splitting angles NFW models are consistent 
with the JVAS/CLASS results in the sense that the NFW predictions are below the upper 
limit put by JVAS/CLASS, for small splitting angles the NFW predictions are far below 
the observational results of JVAS/CLASS. So, both SIS models and NFW models cannot 
explain the observations. We have also tried GNFW models with 1 < a < 2, neither of them 
can explain the observations on both small splitting angles and large splitting angles. All 
the models either predict too many lenses with large splitting angles, or too few lenses with 
small splitting angles. Neither of them can produce the observed large ratio of the number 
of small splittings to the number of large splittings. 

For the same cosmological and lens models we have also calculated the average magni- 
fication bias B = P bs/-P, the results for NFW models and SIS models are shown in Fig. 12. 
We see that, the magnification bias for NFW lenses is bigger than the magnification bias for 
SIS lenses by about one order of magnitude. Though the magnification bias for SIS lenses 
is a constant for all cosmological models, the magnification bias for NFW lenses depends on 
cosmological models and slowly decreases as the splitting angle increases. 

The above results strongly suggest that the halos in the real universe cannot be described 
by a single universal density profile. There must be at least two populations of "halos" in 
the universe: small mass halos with a steep inner density slope, and large mass halos with 
a shallow inner density slope. To test this conjecture, we have calculated the following 
additional lens model: lenses with mass M < M c have the SIS profile (which would produce 
the flat rotation curves seen in galaxies), but lenses with mass M > M c have the NFW profile. 
The critical mass M c for transition from SIS to NFW is determined by fitting our results 
with the JVAS/CLASS observations. Our fitting results indicate that M c ~ lO 13 /i _1 M . 
Interestingly, this critical mass is very close to the cutoff mass of halos below which cooling of 
the corresponding baryonic component will lead to concentration of the baryons to the inner 
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parts of the mass profile (Rees & Ostriker 1977; Blumenthal et al. 1986; Porciani & Madau 
2000). The results for the case with M c = 10 13 h~ 1 M Q are shown in Fig. 13, which fit the 
JVAS/CLASS observations reasonably well. The cosmological models are the same as those 
in Fig. 11. The NFW halos (with M > 10 13 h~ 1 M Q ) have the same concentration parameters 
as those in Fig. 11. Though all the three models are consistent with the JVAS/CLASS 
6" — 15" survey, the SCDM model obviously produces too many lenses with small splitting 
angles. It is interesting that while SCDM produces the most small splittings, it produces 
the fewest large splittings. We have calculated the expected total number of lenses with 
Ad > 0.3" for the JVAS/CLASS sample, the results are 



The expected number of lenses for the LCDM model is remarkably close to the number of 
lenses observed by the JVAS/CLASS survey which is 17 if we exclude the questionable lens 
of B21 14+022 from the sample. 



With the semi-analytical approach we have calculated the probability for forming mul- 
tiple images of a remote source object gravitationally lensed by foreground dark halos. The 
mass density profile of a halo is alternatively taken to be described by the SIS profile, the 
NFW profile, and the GNFW profile. The mass function of halos is assumed to be given 
by the Press-Schechter function. The cosmological model is alternatively the SCDM model, 
the OCDM model, and the LCDM model. Our results show that the lensing probability is 
very sensitive to the density profile of lenses (dark halos). The more the mass distribution 
concentrates toward the center, the higher lensing probability the density profile gives rise, 
which is clearly seen in Fig. 7 - Fig. 10. Compared to SIS lenses, NFW lenses are extremely 
inefficient in producing multiple images. For example, the lensing probability P(> 5") for 
NFW halos is lower than the corresponding probability for SIS halos by more than two orders 
of magnitudes. For GNFW halos, as a decreases, the differences among different cosmolog- 
ical models become prominent (Fig. 8), especially the difference between the SCDM model 
and the OCDM/LCDM model: when a — 1 (i.e. for NFW halos) the lensing probability 
for the SCDM model is lower than the lensing probability for the OCDM/LCDM model by 
more than two orders of magnitudes. This dramatic effect is dominantly caused by the fact 
that halos in a SCDM universe have the smallest concentration parameter (Bartelmann et 
al. 1998) and for smaller a the lensing probability is more sensitive to the concentration 




37, for SCDM 
12, for OCDM 
16, for LCDM 



(89) 



6. Discussion and Conclusions 
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parameter c\. for the NFW case the lensing probability is extremely sensitive to ci, espe- 
cially for small C\ (Fig. 9). Though not so dramatic as the dependence on the inner slope, 
the lensing probability also shows dependence on the outer slope of halos. Comparing the 
right panel of Fig. 8 with Fig. 3, the lensing probability produced by GNFW halos with 
a = 2 is somewhat higher than the lensing probability produced by SIS halos. This mild 
difference is attributed to the different slopes of the a — 2, GNFW profile and the SIS profile 
at large radii (the a — 2, GNFW profile and the SIS profile have the same slope at small 
radii). Similar conclusions are obtained by Wyithe, Turner, & Spergel (2000) and Keeton 
& Madau (2001). Wyithe, Turner, & Spergel (2000) find that the optical depth to multiple 
imaging is a very sensitive function of the profile parameters and the GNFW profile exhibits 
degeneracies between profile parameters with respect to lensing statistics. Keeton & Madau 
(2001) demonstrate that the lensing probability is determined almost entirely by the fraction 
of the halo mass that is contained within a fiducial radius that is ~ 4% of the virial radius. 

Our results also show the dependence of the lensing probability on the cosmological 
parameters, which is particularly manifested in the SIS case. In Fig. 3, for small splitting 
angles the lensing probability in the SCDM universe is higher than the lensing probability 
in the OCDM/LCDM universe, while for large splitting angles the lensing probability in the 
SCDM universe is much lower than the lensing probability in the OCDM/LCDM universe. 
This is explained by the fact that the Press-Schechter function is sensitive to a 8 only for 
large mass halos (Fig. 1). Since small splitting angles are produced by small mass halos, 
and for small mass halos the Press-Schechter function is not sensitive to <r 8 , the lensing 
probability is dominantly determined by the cosmic mean mass density Q m ; thus for small 
splitting angles the lensing probability is highest for SCDM since the SCDM universe has 
the highest mean mass density. Since large splitting angles are produced by large mass halos 
and for large mass halos the Press-Schechter function is (exponentially) sensitive to <r 8 , the 
lensing probability is dominantly determined by <r 8 , thus for large splitting angles the lensing 
probability is lowest for SCDM since the SCDM universe has the smallest <t 8 . For the NFW 
case, in Fig. 6 we see that to produce the same small splitting angle the required NFW halo 
mass is significantly larger than the required SIS halo mass. Thus, for the NFW case <r 8 
takes effect for all A6 > 1", which together with the concentration parameter c\ overtakes 
the effect of Vt m (see Fig. 8). The studies of cluster abundances proposed a correlation 
between Q m and erg as given by equation (22). We have also tested the sensitivity of the 
lensing probability to Q m when equation (22) is satisfied. We see that in general the lensing 
probability varies with Q m ; and for the NFW case the lensing probability is very sensitive 
to a rn (Fig. 10). 

Our numerical results are summarized in Table 1, where we have not included the effect 
of magnification bias (i.e. P is the intrinsic probability), the source object is assumed to be 
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at z s = 1.5, the Hubble constant is assumed to be h = 0.7. To show the sensitivity of the 
lensing probability to the parameters Ug, f2 m , a, and ci, we have evaluated the differentiation 
of lensing probability P(> A9) with respect to cr 8 , Q m , a, and c\ for the LCDM model at 
<7 8 = 0.95, tt m = 0.3, tt A = 0.7, a = 1.5, and Cl (z = 0) = 4 for A# = 5" and A# = 10" 
respectively. The results are 

S ± = 4.9 + 1.7 + 11 

P V ff 8 / n m ,a,ci V ^"1 / a 8 ,a, Cl \ a / a s ,n m , Cl 

+3.1 ( ^) , (90) 

for AO = 5", and 



/ ag,Q m ,a 



% = 5.8 m + 1.8 ("^ + 19 



<Tg,f2 m ,Cl 



+3.4(^1 (91) 



C l / <Tg,Q m ,a 

for = 10". Note the extreme sensitivity of the lensing probability to the slope of the inner 
profile: (dlnP/dlno;) (Tg; Q m)Cl ~ 10. Note also that, if Q m and as enter the lensing probability 
as a combination given by equation (22), then we should expect (din a^/ d\nQ m ) p aci m 7 m 
0.53, but equation (90) and equation (91) show rather smaller values of (din <r 8 / d\nQ m ) p ia ,a 
(0.34 and 0.31 respectively). Thus, in lensing statistics the degeneracy between Q m and a 8 
in equation (22) is broken. This is important since it allows us to separately determine the 
values of Q m and og in principle. 

With the magnification bias being considered, none of the simple models can completely 
explain the JVAS/CLASS observations. The SIS models produce too many large splitting 
lenses, but the JVAS/CLASS observations have detected no lenses with A# > 3" (Helbig 
2000; Phillips et al. 2000). The null result for detecting large splitting lenses are true not 
only for JVAS/CLASS (Phillips et al. 2000), but also for ARCS (Arcminute Radio Cluster- 
lens Search) which is aimed at looking for gravitational lensing events with images separation 
between 15" and 60" (Phillips, Browne, & Wilkinson 2000). While the NFW models produce 
very rare large splitting lenses which is consistent with observations, they produce too few 
small splitting lenses which is against observations since in the JVAS/CLASS survey at least 
17 lenses with 0.3" < A9 < 3" have been discovered (Helbig 2000). None of the models 
can explain the observable large ratio of the number of small splittings to the number of 
large splittings. This strongly suggests that there are at least two populations of halos 
in the universe: small mass halos with a steep inner density slope, and large mass halos 
with a shallow inner slope. We have constructed a very simple two-population halo model 
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to test this conjecture: the mass density of halos with mass < lO 13 /i _1 M is given by 
the SIS profile, while the mass density of halos with mass > 10 13 /i _1 M is given by the 
NFW profile. We find that the results for this model are reasonably consistent with the 
JVAS/CLASS observations. In particular, the number of lenses with A6 > 0.3" predicted 
by the LCDM model is ~ 16, which is remarkably close to the number of lenses observed by 
the JVAS/CLASS survey which is 17. The SIS model predicts too many lenses (~ 37), while 
the OCDM model predicts somewhat too few lenses 12). A similar compound model 
has also been considered by Porciani & Madau (2000) to explain the results from CASTLE 
(CfA-Arizona Space Telescope Lens) survey. 

In summary, with the semi-analytical approach we have shown that the gravitational 
lensing probability is very sensitive to the mass density profile of lenses (especially in the 
central region), the mean mass density in the universe, and the amplitude of primordial fluc- 
tuations. Compared with the observation results of JVAS/CLASS, our calculations indicate 
that the halos in the real universe cannot be described by a single universal universal density 
profile, there are at least two populations of halos in the universe: small mass halos with a 
steep inner density slope and large mass halos with a shallow inner density slope. Ultimately, 
of course, very accurate address of the question is left open by this study. 

We are grateful to E. L. Turner for many helpful discussions and comments, to C. R. 
Keeton and P. Madau for helpful communications. We would like to thank the anonymous 
referee whose comments led to the improvement of this work. This research was supported 
by the NSF grants ASC-9740300 (subaward 766) and AST-9803137. 
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Fig. 1. — Upper panel: The Press-Schechter function / against the velocity dispersion a v of 
dark halos at z = in a SCDM cosmology. The solid curve is for <r 8 = 0.5, the dashed curve 
is for a 8 = 0.6. Lower panel: ( = /JJ^ as a function of a v , for a SCDM cosmology with 
z = and <t 8 = 0.5. Both panels show that for large a v , a poor knowledge in / may give a 
good estimation of os- Since lensing is produced by the high a v part of the distribution, the 
number of lenses sensitively constrains erg. 




z 



Fig. 2. — The differential lensing probability dP/dz for A9 > 5" as a function of lenses' 
redshift z. The source object is assumed to be at z s = 1.5. The lens objects are SIS halos. 
The cosmological models are alternatively: SCDM with Q m = 1 and a"§ = 0.5 (solid line); 
OCDM with Q m = 0.3, fi A = 0, and a 8 = 0.85 (dashed line); LCDM with fi m = 0.3, 
Q\ = 0.7, and as = 0.95 (dotted line). The Hubble constant is h — 0.7. As expected, the 
probable position of the lens is at significantly higher redshifts for the low Q m models than 
for SCDM - note arrows which show median expected redshifts. 
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Fig. 3. — The integral lensing probability P(> A9) as a function of A9. The models are the 
same as those in Fig. 2. 
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Fig. 4. — The lensing equation for the GNFW profile, i.e. equation (62) in the text. Here 
shows the case with a = 1.2 and fi s — 1. The horizontal axis is x, which labels the position 
in the lens plane; the vertical axis is y, which labels the position in the source plane. The 
points (x CT , —y cr ) and (—x cr ,y cr ) satisfy dy/dx = 0. The non-zero roots of y(x) = are ±x . 
Three images are formed when \y\ < y cr , two images are formed when \y\ = y CT , one image is 
formed when \y\ > y CT . So, multiple images are formed when \y\ < y cr , an example is shown 
with the horizontal long dashed line ABC - which has three images: A, B, and C. In the 
paper we consider the splitting angle between the two outside images, i.e the splitting angle 
between A and C. 
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Fig. 5. — S(ci) = l/i](ci) 2 , where c\ is the concentration parameter, r\{c\) = ri/2/r2oo 
is determined by equation (75), where r±/2 is the half-mass radius of a cluster. S(ci) is 
proportional to the surface mass density of a halo at the half- mass radius (eq. [79]). 
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Fig. 6. — Splitting angle A6 as a function of M - the mass of the lens object. The thick 
lines show the splitting angle produced by a NFW lens, the thin lines show the splitting 
angle produced by a SIS lens. The source object is at z = 1.5, the lens object is at z = 0.3. 
The cosmological models are alternatively: SCDM with VL m = 1 (solid lines), OCDM with 
Vt m = 0.3 and fi A = (dashed lines), LCDM with Vt m = 0.3 and fi A = 0.7 (dotted lines). For 
the NFW case, the concentration parameters are alternatively 5 for SCDM, 9 for OCDM, 
and 7 for LCDM. Notice the enormous sensitivity at small splittings and small (e.g. galactic) 
mass halos to the steepness of the inner part of the mass profile. 
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Fig. 7. — The differential lensing probability dP/dz for A6 > 5" as a function of lenses' 
redshift z. The source object is at z s = 1.5. The lens objects are alternatively NFW halos 
(left panel, i.e. GNFW halos with a = 1), GNFW halos with a = 1.5 (central panel), 
and GNFW halos with a = 2 (right panel). The concentration parameters are alternatively 
given by equation (76), equation (77), and equation (78). The cosmological models are 
alternatively: SCDM with Q m = 1 and cx 8 = 0.5; OCDM with Q m = 0.3, Q A = 0, and 
a 8 = 0.85; OCDM with Vt m = 0.3, fi A = 0.7, and a 8 = 0.95. The Hubble constant is h = 0.7. 
The extreme sensitivity of lensing to the slope of the inner profile is evident. Arrows show 
median expected redshifts. 
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Fig. 8. — The integral lensing probability P(> Ad) as a function of A9. The models are 
the same as those in Fig. 7. Note the relatively slow fall-off of probability with increasing 
splitting angle for the OCDM and LCDM cosmologies: for all a the number of 10" splittings 
expected is not very much smaller than the number of 1" splittings expected. 
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Fig. 9. — The integral lensing probability P(A6 > 5") as a function of lenses' concentration 
parameter c\. The cosmological models are the same as those in Fig. 7. The source object 
is at z s = 1.5. The lens objects are alternatively NFW halos (left panel), GNFW halos with 
a = 1.5 (central panel), and GNFW halos with a = 2 (right panel). 
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Fig. 10. — The integral lensing probability P(A0 > 5") as a function of Q m when Q m and 
as are constrained by equation (22). The lens models are the same as those in Fig. 7. The 
source object is at z s = 1.5. The cosmological models are OCDM (solid lines) and LCDM 
(dashed lines). The Hubble constant is h — 0.7. Note that, for models normalized to give 
the correct z = cluster abundances, the sensitivity of the lensing probability to the halo 
profile shape is far stronger than to the matter density tt m . 
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Fig. 11. — Comparison of our semi-analytical results with the JVAS/CLASS survey. The 
JVAS/CLASS results are shown with the histogram. The null result for lenses with 6" < 
A6 < 15" is shown with the horizontal line with a downward arrow indicating that is an upper 
limit. Our semi-analytical results (allowing for amplification bias) are shown with solid curves 
(SCDM with tt m = 1 and a 8 = 0.5), dashed curves (OCDM with VL rn = 0.3, fi A = 0, and 
(7 8 = 0.85), and dotted curves (LCDM with Q m = 0.3, £1\ = 0.7, and <r 8 = 0.95). Two groups 
of lenses models are shown: SIS lenses, and NFW lenses with the concentration parameter 
ci = 5 for SCDM, 9 for OCDM, 7 for LCDM. The NFW+SCDM lensing probability has 
been multiplied by a factor 10 L5 to fit it on the same scale as other models. For comparison 
we have also shown the result for a NFW+LWDM model with the dashed-dotted curve, 
which has the same parameters as the NFW+LCDM model except that C\ = 3.5. 
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Fig. 12. — The average magnification bias as a function of splitting angles. The solid, 
dashed, and dotted lines are for NFW lenses, the thin dashed-dotted line is for SIS lenses. 
The cosmological and lens models are the same as those in Fig. 11. The magnification bias 
for SIS lenses is a constant independent of cosmological models. 
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Fig. 13. — Fitting of the two-population models to the JVAS/CLASS results. The cos- 
mological models are the same as those in Fig. 11. The lens models are: SIS halos with 
M < 1O 13 /i _1 M , NFW halos with M > 10 13 /i _1 M . The concentration parameters for 
NFW halos are alternatively ci = 5 for SCDM, 9 for OCDM, and 7 for LCDM. 
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Table 1. Summary of the Results for the Lensing Probability 



Models 


lgP(2") 


lgP(5") 


lgP(10") 


lgP(20") 


lgP(40") 


SCDM a (a = 1) 


-8.16 


-8.54 


-9.17 


— 10.4 


— 12.8 


SCDM a (a = 1.5; 


) -5.14 


-5.46 


-5.98 


-6.96 


-8.84 


SCDM a (a = 2) 


-2.82 


-3.10 


-3.49 


-4.17 


-5.36 


SCDM a (SIS) 


-3.00 


-3.33 


-3.78 


-4.56 


-5.97 


OCDM b (a = 1) 


-5.64 


-5.68 


-5.79 


-6.04 


-6.60 


OCDM b (a = 1.5) -4.37 


-4.49 


-4.66 


-5.00 


-5.62 


OCDM b (a = 2) 


-2.81 


-2.98 


-3.20 


-3.56 


-4.17 


OCDM b (SIS) 


-3.23 


-3.45 


-3.74 


-4.23 


-5.04 


LCDM c (a = 1) 


-6.12 


-6.17 


-6.28 


-6.54 


-7.09 


LCDM c (a = 1.5) -4.50 


-4.62 


-4.80 


-5.14 


-5.77 


LCDM c (a = 2) 


-2.74 


-2.90 


-3.12 


-3.48 


-4.08 


LCDM C (SIS) 


-3.04 


-3.24 


-3.51 


-3.96 


-4.69 


LWDM d (ct = 1) 


-7.45 


-7.46 


-7.50 


-7.59 


-7.80 


a SCDM: n m = 


1, Qa = 0, CT 8 


= 0.5, ci( 


z — 0) = 5 for a — 1. 




b OCDM: tt m = 


0.3, Q A = 0,a 8 = 0.85, 


Cl (z = 0) = 


9 for a — 1. 




C LCDM: VL m = 


0.3, Q A = 0.7, 


a 8 = 0.95, a(z = 0) : 


= 7 for a = 1. 




d LWDM: Q m = 


: 0.3, tt A = 0.7, a 8 = 0.95, a(z = 0) 


= 3.5 for a = 


: I. 



Note. — The GNFW parameter a is defined in equation (60). SIS = singular 
isothermal sphere. The source object is at z s = 1.5. The Hubble constant is 
h = 0.7. The magnification bias is not included, so P(A9) = P(> A9) is the 
intrinsic lensing probability. 



